Laws of Large Numbers for the Occupation Time of an 
Age-Dependent Critical Binary Branching System 



J. Alfredo Lopez-Mimbela 
j alf redoQcimat . mx 

Antonio Murillo-Salas 
murilloOcimat .mx 

Centro de Investigacion en Matematicas, Guanajuato, Mexico 

Abstract 

The occupation time of an age-dependent branching particle system in R'' is considered, 
where the initial population is a Poisson random field and the particles are subject to symmetric 
a-stable migration, critical binary branching and random lifetimes. Two regimes of lifetime 
distributions are considered: lifetimes with finite mean and lifetimes belonging to the normal 
domain of attraction of a 7-stable law, 7 G (0, 1). It is shown that in dimensions d > a-f for the 
heavy-tailed lifetimes case, and d > a for finite mean lifetimes, the occupation time proccess 
satisfies a strong law of large numbers. 
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1 Introduction and background 

In this paper, we obtain laws of large numbers for the occupation time process of a random popula- 
tion living in the d-dimensional Euclidean space M''. The evolution of the population is as follows. 
Any given individual independently develops a spherically symmetric a-stable process during its 
lifetime r, where < a < 2 and r is a random variable having a non-arithmetic distribution 
function, and at the end of its life it either disappears, or is replaced at the site where it died by 
two newborns, each event occurring with probability 1/2. The population starts off from a Poisson 
random field having Lebesgue measure A as its intensity. We postulate the usual independence 
assumptions in branching systems. 

Two regimes for the distribution of r are considered: either r has finite mean /i > 0, or r 
possesses a distribution function F such that F(0) = 0, F{x) < 1 for all x E [0,oo), and (with 
g{u) ~ h{u), as u ^ 00, meaning g{u)/h{u) — > const, as u ^ 00) 

:= 1 -F(u) ^ n-Tr(l -7)"^ as u — > 00, (1) 
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where 7 £ (0,1) and F denotes the Gamma function, i.e., F belongs to the normal domain of 
attraction of a 7-stable law. In particular, this allows to consider lifetimes with infinite mean. 

Let X{t) denote the simple counting measure on whose atoms are the positions of particles 
alive at time t, and let X = {X{t), t > 0}. When r has an exponential distribution it is well 
known that the measure-valued process X is Markov. In the literature there is a lot of work about 
the Markovian model. Our objective here is to investigate the case when r is not necessarily an 
exponential random variable, in which case {X(t), t > 0} is no longer a Markov process. Another 
striking difference with respect to the case of exponential lifetimes arises when the particle lifetime 
distribution satisfies ([T]). When the distribution of r possesses heavy tails, a kind of compensation 
occurs between longevity of individuals and clumping of the population: heavy-tailed lifetimes 
enhance the mobility of individuals, favouring in this way the spreading out of particles, and thus 
counteracting the clumping of the population. Since clumping goes along with local extinction (due 
to critical branching), a smaller exponent 7 suits better for stability of the population. As a matter 
of fact, Vatutin and Wakolbinger [IB] and Fleischmann, Vatutin and Wakolbinger [S] proved that 
X admits a nontrivial equilibrium distribution if and only if d > 7a. This contrasts with the case 
of finite- mean (or exponentially distributed) lifetimes, where the necessary and sufficient condition 
for stability is d > a. As we shall see, such qualitative departure from the Markovian model 
propagates also to other aspects of the branching particle system, such as the large-time behavior 
of its occupation time. 

Recall that the occupation time of the measure-valued process X is again a measure-valued 
process J = {Jt, t > 0}, which is defined by 



theorems for occupation times of branching systems have been extensively studied in the context of 
exponentially distributed lifetimes. Cox and Griffeath [B] proved a strong law of large numbers for 
the occupation time of a critical binary branching system. Moreover, in [6j it is proved a central 
limit-type theorem for the occupation time of the critical binary branching Brownian motion. 
Meleard and Roelly [13] extended the law of large numbers of [6] to branching populations with 
general finite-variance critical branching, and quasi-stable particle motions. Bojdecki, Gorostiza 
and Talarczyk [3 IH [S] have investigated the limit fluctuations of the re-scaled occupation time 
{Jxit) := JiT^t > 0} of branching systems, T being a parameter which tends to infinity. They 
have shown that these limits are processes which exhibit long-range dependence behavior, such as 
fractional Brownian motion and sub-fractional Brownian motion. See also [1] for related results, 
where the underlying process is a branching random walk in the d-dimensional lattice. 

In this paper we will prove that, in dimensions d > 07 for heavy-tailed lifetimes, and d > a for 
finite-mean lifetimes, the occupation time of the process X satisfies a strong law of large numbers. 




for all bounded measurable function tp : 



M+, where the notation {ip^v) means J ipdv. Limit 
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Namely, a.s. for any positive continuous function ip with compact support, 

t~^{'il^,Jt) — > {ip,A) as t — > oo. 

Also, we prove that in dimensions d < 07 for heavy-tailed lifetimes, and d < a for finite-mean 
lifetimes, the normalized occupation time t~^Jt converges to zero a.s. in the sense that, with 
probability 1, for any ball A C M'^ of finite radius, 

'^{Xs(A)>o} ds — >0 as t — ^00. (2) 

Jo 

These results complement — and partially extend — those of [B] and [13]. We point out that 
dimension-dependent behaviors, or parameters, are a typical characteristic in this theory because 
properties of the branching system treated here are highly related to the transience-recurrence be- 
havior of the particle motions. Notice also that, in contrast with the case of finite-mean lifetimes, 
in the presence of heavy-tailed lifetimes the dimension restriction varies according to the decay 
exponent of the tail. This phenomenon is reminiscent of the interplay of population clustering and 
longevity of individuals quoted above. 

Our proofs use techniques from [B], [U] and [13]. To prove the strong law of large numbers in 
case of heavy-tailed lifetimes, we first consider the case of "intermediate dimensions" 07 < d < 2a, 
and deal afterward with "large dimensions" d > 2a. Aiming at applying the Borel-Cantelli lemma, 
in case of intermediate dimensions we use the re-scaled occupation time process to upper-bound the 
variance functional of the occupation time. This step employs certain Fourier-transform techniques 
that we adapted from [3]. We were unable to extend this method to dimensions d > 2a due to 
the lack of proper upper-bounds for the variance functional of the re-scaled occupation time. To 
deal with the case of large dimensions we follow the approach of [13]. We use a Markovianized 
branching system, introduced in Section 4.1 below, which allows us to directly apply the well- 
known self-similarity of the symmetric a-stable transition densities. We remark that, in order to 
use this procedure, we need to assume that the lifetime distribution possesses a continuous density 
function. This contrast with the case of low dimensions, where no absolute continuity condition is 
required. We think, however, that the result should be true for a general lifetime distributions. 

In case of a general non-arithmetic lifetime distribution having finite mean, our proof of the law 
of large numbers is carried out using estimates for the variance functional of the occupation time 
process, as well as bounds for the a-stable transition densities. The almost sure convergence ([2]) is 
proved by combining Borel-Cantelli's Lemma with some estimates from |16j related to extinction 
probabilities. 

The analysis at the "critical dimensions" d = 07 in the heavy-tailed case, and d = a for finite 
mean lifetimes, is much more difficult to carry out, as can be seen from where the occupa- 
tion time (at the critical dimension d = a/ 13) of the so-called (d, a, /3)-superprocess is considered. 
The approach there strongly relies on the classical semilinear equation characterizing the Laplace 
functional of the occupation time, see Lemma 3.4 in In our case, due to the non-exponential 
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lifetimes, we do not have the above-mentioned equation. Thus, laws of large numbers for our model 
at critical dimensions remain to be investigated. 

2 Laws of large numbers 

Following [3], we define the re-scaled occupation time process {Jrit) ■= JtT, * > 0}, i.e., for any 
positive bounded measurable function ip, 

{if, Jrit)) := [ {^,Xs)ds = T [ {ip,XsT)ds, t>0. (3) 
Jo Jo 

Notice that, by Fubini's theorem, 

E{^,Jt{1)) = {v,A)T, (4) 

since IE((/7, Xt) = {ip, A). We remark that studying the asymptotic behavior of ((/?, Jt) /t as t — > oo, 
is the same as investigating the asymptotic behavior of {ip, Jt{^))/T as T — > oo. 

In what follows, C^(M'^) denotes the space of nonnegative continuous functions ip : M*^ — > M+ 
with compact support. The main results of this paper are the following theorems. 

Theorem 2.1 Let F be a non- arithmetic distribution function satisfying ([ip. 

(a) Assume that aj < d < 2a. Then, a.s. for any ip G C^{M.'^), 

T'\ip,JTil)) — >{'P,A) as T — ^oo. (5) 

(b) Suppose that d > 2a and F possesses a continuous density f. Then, a.s. for all ip & 

T-\ip,JTil)) — >{'P,A) as T — ^oo. (6) 

Our next theorem complements the law of large numbers of [6] and [13], which were proved only 
in the case of exponentially distributed lifetimes. 

Theorem 2.2 Assume that d > a, and let F be a non- arithmetic distribution function with finite 
mean /j, > 0. Then, a.s. for any ip € C^(M'^), 

T-\ip,JT{l)) ^ {ip,A) as T^oo. (7) 

Theorem 2.3 Let F be a non- arithmetic distribution function which satisfies (1^, and assume that 
d < aj. Then, a.s. for any ball A CM.'^ of finite radius, 

'^{Xs{A)>o}ds — >0 as T — > oo. (8) 

Jo 
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Remark 2.4 (a) Notice that condition 07 < d < 2a in Theorem \2.1\ allows d < a, which contrasts 
with the classical case of exponentially distributed lifetimes, where d > a. 

(h) When the particle lifetimes have an exponential distribution with mean , Theorem \ 2. 2\ reduces 
to Theorem 4 of \13^ . 

(c) In case of low dimensions d < 07, a genuine counterpart to Theorem \2. 1\ would be a statement 
ensuring a.s. vague convergence of Jt{1) to the zero measure as T ^ 00. This was proved by 
Sawyer and Fleischman ]15^ for the occupation time of critical branching Brownian motion (see also 

for a related result regarding super- Brownian motion's occupation time). For our model, here 
we prove only the slightly weaker result which implies that, with probability 1, the proportion of 
time that the branching system charges any given bounded set vanishes asymptotically as T —> 00. 
The specific form ^ of a.s. convergence was suggested to us by an anonymous referee. 

(d) The extent of dimensions in our results is narrow due to our choice of critical, binary-branching 
mechanism. A less restrictive assumption, such as critical (1 + fi) -branching, (5 € (0,1), would 
expand the dimension range. 



3 Some moment calculations 



Let Zt{A) denote the number of individuals living in ^ € B{M.'^) at time t, in a population starting 
with one particle at time t = 0. Following [TT] we define 



Qtfix) := 



1 - e-<^'^'> 



(9) 



where ip G C+(R'^) and E X means that the initial particle is located at x G M*^. Since the initial 
population Xq is Poissonian, we have 



]Eg-('^,Xt> ^ gxp ( - / E^. [l - e-^'^'^'H dx 



exp 



QMx)dx] , if G C+(M°'). 



(10) 



Let {rfc, /c > 1} be a sequence of i.i.d. random variables with common distribution function F, and 
let 



k=l 



where the random sequence {Sk, A; > 0} is recursively defined by 

So = 0, Sk+i = Sk + Tk, k> 0. 
For any p = 1,2, < tp < ip_i, . . . , ii < 00, ipi,ip2, . . . ,ipp G C^{M.'^) and 0} 



, . . . ,wp 



G M, we 



define t = (ti,t2, • • • ,tp), t - s = {ti - s,t2 - s, . . . ,tp 



Hp) 



9i,...,9py and 



Qf0(p)(x)=E, 
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Let {Bg, s > 0} denote the spherically symmetric a-stable process in M*^, with transition density 
functions {pt{x,y), t > 0, x,y € M"'}, and semigroup {St, t > 0}. Our moment calculations use the 
following result which is borrowed from [TT| (Section 4.3), and which we include for convenience. 



Proposition 3.1 The function Q-^ 6 (^p^ satisfies 



p-i 



dN, 



As in (jlOp . since the initial population is Poissonian we have 



E 



exp 



exp 



E. 



dx 



Q^9(^p){x)dx 



(11) 



Using criticality of the branching, and that Lebesgue measure is invariant for the semigroup of the 
symmetric a-stable process, it is easy to see that 

m{t,ip) ■.= E[{^,Xt)] = {ip,A), t>0, ip£C+(R'^). 



(12) 



Lemma 3.2 Let < s < t < oo and ip^ip £ C+(IR'^). Then, 
C^{s,ip;t,tp) := E^. 

= E, 

where mx{t,Lp) = E^.[(99, Zf)] . 



Lp{Bs)'4){Bt) + / ruB^it - r,i;)mB,{s - r,ip)dNr 




(13) 



Proof: In order to preserve the notation in Proposition 13.11 we put p = 2, ti = t, t2 = s, (fi = ip 
and ip2 = (f. Then we have 

92 



ei =6(2=0+ 



where 



Qlei2){x) = E^ 



't2 Q 



' -§^QtrOi2)iBr)^QtrSi2){Br)dNr 

QtA2){Br)dNr 



^2 



^ iQtA.^m) ^^^^^^ 



-^2{Bt,)e-'^'(^^2'^ J^'' {Ql_MBr)) -^Ql_MBr)dNr 
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Evaluating at = ^2 = we finish the proof. 

Proposition 3.3 Let < s < t < oo and '4>,ip £ C+(E'^). Then, 



□ 



Jo 

where U{r) = E^o^*''(0- 

Proof: We put p = 2 in (jlip and use the same notations as in the proof of Lemma |3.2[ Then, 

92 



E[{^i,Xt,){ip2,Xt,)] 



ddidd2 



exp 



QlO[2){x) dx 



ei=e2=o+ 



Qf6(2)ix)dx 



861862 

+ I -§Q-Ql^{2){x) dx I Aq?0(2)(x) dx 



61=62=0+ 

Cx{ti,fi;t2,(f2)dx+ I m.x{ti,(fi)dx I mr,{t2, (p2)dx, 



and from Lemma 13.21 we obtain 

Cis,ip;t,ij)= [ 
which completes the proof. 



^p{Bs)ip{Bt) + / niB.it - r,ij)mB,{s - r,ip)dNr 
■Jo 



dx. 



(15) 
□ 



4 Markovianizing an age-dependent branching system 

In this section we introduce a Markovian branching system {Xt, t >0} which will be used to prove 
Theorem l2. 11 (b). Let X = {X^, t > 0} be the branching system defined in Section 1. For any t > 0, 
let Xt denote the population in x M"^ (M+ = [0,oo)) obtained by attaching to each individual 
6x € Xt its age. Namely, for each t > 0, 



(16) 



where rjl and Q denotes respectively, the age and position of the i^^ particle at time t, and the 
summation is over all particles alive at time t. Let us assume that Xq is a Poisson random field 
on M-|_ X M'^ with intensity measure F x A. Here, F also means the Lebesgue-Stieltjes measure 
corresponding to F. The probability generating function of the branching law is denoted by 
Thus, for critical binary branching, $(s) = ^(1 + s^), —1 < s < 1. 



Given a counting measure i' on 
define 



and a measurable function 



(0,1], we 



G4,{iy) := exp ((log </>, z^)) . 
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It can be shown that the infinitesimal generator of {X^, t > 0} evaluated at the function G^{i') is 
given by 

eG,M=G4.)(£*M±^<«^>^^..). (17) 

where 

X{u) = -I^ u>0, (18) 
1 — r [u) 

is the hazard rate function associated to F, and 

C(l){u, x) = ^^^^ + Aa^{u, x) - A(n) [0(n, x) - 0(0, x)] , (19) 

where the function (f> is such that </>(•, x) G C^(M+) for any x G W^, and (j){u, •) G C^{W^) for any 
G M+. Here C^(M+) denotes the set of all bounded functions with continuous first derivative, 
and C^{R'^) denotes the space of infinitely differentiable functions from M"^ to M, having compact 
support. The operator £ is the infinitesimal generator of a Markov process on x R'^ whose 
semigroup is denoted by {Ti,t > 0}, see [H] for details. 

Proposition 4.1 Let X = {Xt, t > 0} as before and let Xq be a Poisson random field on x M*^ 
with intensity measure F x A. The joint Laplace functional of the branching particle system X and 
its occupation time is given by 



E 



for all measurable functions ip, (f> : M+ x M'^ — > with compact support, where V^^cp satisfies, in 
the mild sense, the non-linear evolution equation 

^vU{n,x) = CVU{u,x)-\{umi-VU{^,^))-{-^-VU{^,m 

+<l){u,x){l-VU{u,x)), (20) 
V^(t>{u,x) = i-e-'^^"'^). 

Proof: The proof is carried out using the martingale problem for {Xt, t > 0}, and Ito's formula. 
We omit the details. 



5 Proof of Theorem 12.11 

We shall prove the law of large numbers in two steps. First we show that the result holds for 
intermediate dimensions 07 < d < 2a; this part of the proof relies on the non-Mar kovian branching 
system defined in Section 1, and uses upper bounds for the covariance functional. In the second 
step, we consider "large" dimensions d > 2a, and in this case we use the Markovianized branching 
system described in Section HI 
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5.1 Proof of Theorem [M] (a) 

In this section we assume that 07 < d < 2a. 

Lemma 5.1 Suppose that the hypothesis in Theorem \2. 1\ hold. Then, for each e > and allT > 
large enough, 

P {\T~\ip, Jt(1)) - (95, A)| > e) < I (c3r-2 + ciT-i + csT-^/" + C4T^-'^/°) , 
for some positive constants ci , . . . , C4 . 

Proof: Let e > be given. Then, using Chebyshev's inequahty and 

P(|r-i(^,JT(l))-(^,A)| >e) < ^E{T-\^,Ml))-{ip,A)f 

rCov((^,JT(l)),(^,JT(l))) 



1 



e 
1 



1 /■! 

Gov (((/?, X^t), {(p,XtT)) dt ds, 



'0 ^0 

where the last equahty foUows from By changing the order of integration we obtain that 

P{\T-\ip,JT{l))-{^,A)\>e)<\ C dv f CoY{{^,XuT),{v,X,T))du. (21) 

e Jo Jo 

Therefore, from Proposition 13.31 we deduce that 

P{\T'\ip,JT{l))-{v,A)\>e) < (/) + (//), (22) 

with 



and 



e Jo Jo 

{II) := — dv du I dU(Tr){ipST{v+u-2r)^,A), 
e JO JO JO 

where, to obtain (II), we used self-adjointness of St with respect to A, t > 0. Our next goal is to 
derive useful upper bounds for the two integrals (I) and (II). Firstly, by performing the change of 
variables s = {v — u)T and t = vT, we get that 



-(/) = ^ dt ds{ipS,ip,A) 



dt I ds{ipSs^p, A) 
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for any ^ > 0, where 



{ipSsip,A)ds =111 '^{x)ps{x - y)ip{y)dy dxds 







(f{x)(p{y) ps{x - y)ds dy dx 
Jo 

< I I ^ix)ip{y)const.(\x-y\'^-'^ + t^-'^/'')dydx 



since 



^ Ps{x - y)ds < const. (\x - yl""'^ + t^"'^/" 
due to self-similarity of the a-stable transition densities. Notice that 

ip{x)Lp{y)\x — y\°'~'^dydx < oo, 



which, for d > a, follows from Lemma 5.3 in [9]. Hence, 

[ dt [ {ipSsip,A)ds < const. / / ip{x)ip{y)[ (\x-y\"'-'^ + t^''^^°')dtdydx 
J A Jo Jr'^ Jr'^ Ja^ ' 

= ci(T- A) + c2(ri-^/° -^1-'^/") 

for some constants ci,C2 > 0. Therefore, 



lie T j^2—d/a \ 

(/)<3UI + ci^ + C2^^ , (23) 



where 



C3 = dt {ifSs^p, A)ds. 
Jo Jo 



Before estimating the integral (//), we recall that U{t) ~ t'^/T{l + 7) as t — > 00 because of 
F{t) ~ t~'^/r(l — 7), see [2], p. 361. Then, writing ip for the Fourier transform of 93, we obtain 



g2 /■! rv j-u 



JO 



{II) = dv du dU{Tr)—-. / d2/|(^(2/)|2e-^('^+"-2'-)l?^ 



f\v [du i/yW{y)? [ dU{Tr)e-^^-+--^'^\y\' 
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dv du dy\ip{y)\'^ / e-^(^+"-2'-)l2^l\T-idr, 

Jo iR<* Jo 



r(l + 7)(27r)rf 

and, after the change of variables z = {T{v + u — 2r))^/"y, we conclude that 

-{II) ~ , ^^J, C dv r du [ dz r T-'^/^{v + n - 2r)-'^/" 
2^ ' T{l + j){27Ty Jo Jo Jr'^ Jo ^ ' 

X |(^(T-'^/"('t; + u- 2r)-'^/"z)|2e-l^l"r^-i(ir 

< / dze-\^\- r riv + u-2r)-''/^r^-'dr, 

- r(l + 7)(27r)^ A. Jo Jo Jo ^ ' 
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where to obtain the last inequality we have used the well known fact that < {2tt) '^{\ip\,A) 

for any L^-function Lf. Changing the order of integration into the above expression yields 

r(l + 7)(2^)'^ 111^ ' 



Jo 



dze ' ' / av r ' ar 



r(l + 7)(27r)'^ J^, Jo Jo 1-d/a 

r(l + 7)(2^)'^ 1-d/a J^, Jo Jo 

Notice that in the above calculations we have imphcitly assumed that a. The case d = a can 
be treated in a similar way (and renders the same conclusion (|24p). Changing again the order of 
integration we get 



Jr 



r(l +7)(27r)'^ 1-d/a 



r(l + 7)(2,r)<' (1 - (i/o)(2 - rf/Q) 7„ 
where the last equality is finite since by assumption d < 2a. Hence, for T large enough 

< ^C4r^-'^/". (24) 

Therefore, 

P {\T-\^, Jr(l)) - (95, A)| > e) < I (cgT-^ + ciT~' + caT"'^/" + c^T''-^/^') . 



□ 



Proof of Theorem 12.11 (a): Let e > and a > 1 be given constants, and let T„ = a" for 
n = 1, 2, . . . . Then, 



£p(|r„i(^,Jr„(l))-(^,A)|>e) 



n=l 

00 

00 



2 



n=l 



due to the assumption d > 7a. Therefore, for any given tp € C^(M'^), a.s., 

^ (V'; '^r„ (1)) — > if, A) as n — ^ 00. 
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Now we observe that, for each T > 1, there exists some non-negative integer n{T) such that 
a"(T) < 7^ < a"('^)+i, and that n{T) — ^ oo as T — > oo. Hence, 

{V,JaMT)il)) (y,^T(l)) ((/?, J^n(T) + l(l)) 

aniT)+l - J, - ^n(T) 

and 

< limmr — < limsup — < ((/?, A)a, 

these inequahties being true for any a > 1. Letting a ^ 1 yields that 

hm r-i(^,Jt(l)) = (<^,A) 

T — ^oo 

a.s., where the null set (where the hmit may not exist) depends upon (p. Nonetheless, a nuh set 
can be chosen not to depend on (p as is the proof of Theorem 1 in [lOj. □ 

5.2 Proof of Theorem EH] (b) 

Throughout this section we assume that d > 2a. The proof of part (b) in Theorem 12.11 fohows. as 
in part (a), from Chebyshev's inequality 

p{ '^^'^^^;^^'^^' >e}<^Var(^,J,), t>0, e>0 (25) 

and Lemma 15.21 below. Recall that A is defined in (1181). 



Lemma 5.2 i) Let cf) -.W^ ^ M_|_ he a measurable function with compact support. Then, for each 
t > 0, 

E(<^,J4) = A)t, (26) 

and 

Var((/), Jt) < (A,F)Const(0)(t + i^-d/") + 2Const((/))(i + t^-'^/"). (27) 

Proof: First we prove (j26p . For the given function (j) we define the extended function (j){u,x) = 
(/)(x), (n, x) G M-i- X W^. Then, for any A; > we define 



^-{Vtik^),FxA) 



where Vt{k(j)) satisfies (|2U|) with substituted by k(j), and -0 = 0. Notice that 

E{cP,Jt) = -A.]E[eKp-{k^,Jt)]\k=o+ 

^V^\k^), Fx a) exp (- {V,'{k4>), Fx A)) 

Thus, putting Vt4> := jrV^{k^) 1^=0+ and recalling that V^°(0(^) = 0, we obtain 
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Vo(j){u,x) 



£Vt(j){u,x) + (t){u,x) 
0, 



or 



Vt(t){u,x)= / ft^s4>{u,x)ds = / St-s(l){x)ds. 
Jo Jo 

Consequently, using that A is invariant for the a-stable semigroup, 



(28) 



This proves (I26p . The proof of (I27p goes as fohows. Define (p as before. Differentiating Vt(A;</>) with 
respect to k and using equation (pOj) . we obtain 



^Vt{k^){u,x) = C^Vt{k^){u,x)+^{u,x){l - Vt{k$){0,x)) - k^{u,x)^Vt{k^){u,x) 



-AN 



-<I>'(1 - Vtik^iO,x))—Vtik^){0,x) + — y,(A:(^)(0,x) 



and 



^3 - 52 - _ a - _ 92 - 

■^^^Vt{k(l)){u,x) = C^Vt{k(t)){u,x) - 2(j){u, x)—Vt{kcl)){u, x) - k(l){u, x)-^Vt{k(p){u, x) 



-A(n) 



^" {1-Vtik^)i0,x)) [—Vtik^)iu,x 



d 



-$'(1 - Vt{k^){0,x))-^Vtik^){0,x) + ^y,(fc<^)(0,x) 

Letting l^c/) = ^Vj(A:(/))|fc=o+ we get that 

^VtHu,x) = CVt4>{u,x) - X{u)^"{l) (vMO,x)y - 2^{u,x)Vt^{u,x) 



(29) 



From ([28]) and we obtain 
Vi<^(u,x) 

Note that Var((/>, Jt) = -{Vtct){*, •), F x A). Therefore 











/o 





Var((/), Jt] 



A(*) / f;0(*,.)(ir +20(*,.) / T;(/.(*,.)dr 



ds,F X A 



A(*) ( / 7;0(*,.)dr) ,FxA)ds 



+2 ('^^*'*^/ f;0(*,«)dr,F X A^ (is 
(^) + (i?). 



(30) 
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Notice that, under the choice of Tt(f){u,x) = St(p{x) for all t>0, and that (A,F) < oo. In fact, 
using that A(n) ~ u^^ and f(u) ^ u~'^~^, we get that, for A > sufficiently large. 



oo 



(A,F) = / \{u)f{u)du 





A POD 

X{u)f{u)du+ / X{u)f{u)du 
J A 

A POO 
X{u)f{u)du + / u^^u^^^^du 
J A 



< oo. 

Now, 

i^) = J^ (^^*^ (i) ^^^^'^^'^^ ,F X ds = {X,F) (^(^J^ SrHr^ ds. 
Also, it can be shown that 

^ ^(^^ Srcpdr^ < Const(02)(t + t^~'^/"), 

and consequently, 

(A) < (A,F)Const(</.2)(t + t^"'^/"). 



Similarly, for the second term in (|30 



ft 

{B) = 2l (</>(*,•)/ Tr(l){*,»)dr,F X A) ds = 2 I (/ Sr())dr,A)ds, 



\J0 



where 

ft / /.s 



^ Sr(l)dr,A^ ds < Const(</.)(t + ^2-'^/"), 



hence, 

(S) < 2Const(0)(t + t2-'='/"). 
Finally, combining the bounds for (A) and (B) we obtain the result. □ 

6 Proof of Theorem D 

Suppose that is a non-arithmetic distribution function supported on the non-negative real line 
and having finite mean > 0, and let d > a. As in the proof of Lemma |5. 11 we have that 

P{\T-^{ip,JT{l))-{^,A)\>e) <^ C dv r Coy {{if, XuT),{v,X,T))du. (31) 

e Jo Jo 
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Therefore, due to Proposition 13.31 

P{\T-\^,Jt{1)) - (v^,A)| > e) < (/) + (//), 

where 

W •= 4 / / du{ipST{v-u)f,-^), 

e JO Jo 



(32) 



and 



(1/):=-/ / / ((5T„-r9')(5T^,-r<^),A)c^[/(r)du(^^;. 
e JO JO JO 

We recall the bound ()23p for (/). It remains to upper-bound (//). Performing the change of 
variables h = r/T in (//), and using the elementary renewal theorem (see e.g. [12j, p. 188), we 
have that, for T large enough. 



ill) 



1 fV fU 



JO JO 



{{ST{u-h)V){ST(v-h)V), 



U{Th) 
Th 



Th 



du dv 



1 fV fU 



2T 

e^/^ Jo Jo Jo 

/■! rii rv 



{ {STiu-h) ^) {ST{v-h) ip),A)dhdu dv 



{{ST{u-h)^){ST{v-h)^),^)dudhdv. 

e ^ JO JO Jh 

After performing several changes of variables one can see that, for all T large enough. 



ill) 



< 



Jo 

2 



V ft 



{Ssip){x){Stip){x)ds dt dx dv 



T 



2 



Jo 



JO 



{Ss'^){x){St^){x)ds dt dx dv 



T 



On the other hand, one can show, as in ^13j, that 



Jo 



(p{y)ip{z) / pt+s{y - z)dsdtdydzdv. 



1^ j'^ pt+s{y - z)ds dt<c{\y- zp-^ + v^-'"^) 
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for some constant c > 0. Hence, for any fixed ^ > 0, and all T large enough, 



(11) < 



+ 



+ 



ip{y)ip{z) / / pt+s{y - z)ds dt dy dz dv 
Jo Jo 

ip{yMz)\y - z\^''-''dydz{T - A) 



< 



ip{y)ip{z)dy dz 
f{y)<f{z 



V fV 



+ 



e2^r 



H — T^c 



3 — d/a 
Pt+s{y — z)ds dt dy dz dv 

ip{y)ip{z)\y - zl^^'^'^dydz 



JO 



(f{y)(f{z)dydz- 



3 — d/a 



The proof concludes with an application of Borel-Cantelli's Lemma, using that d/a > 1, (j32p . and 
the bounds for (I) and (II). 

7 Proof of Theorem ESI 

In this section we assume that d < 07, and that F is a distribution function satisfying ([1]). Arguing 
similarly as at the end of the proof of Theorem 12. 11 Lemma |7. II below yields the theorem. 

Lemma 7.1 Let A C^"^ be a ball. Then, for all e > 0, and for all t sufficiently large, 

P {t-' 1^ lix.iA)>o}ds > < (1 - e-)-^ ^,t-W»^^y'Hi^'^'^/» + ^it"^) , (33) 
for some positive constants c and Ci . 



Proof: Notice that, by Markov's inequality, 

P (^t-^ ^ l{X4A)>o}ds > < (1 - e-T^E 
Moreover, since the initial population is Poissonian, 



1 - exp <^ t^^ / 1 



^{Xs{A)>0} 



ds 



(34) 



E 



* ^ Jo ^{Xs(A)>0}d.S 



E 



exp 



1 _ q-Jq l{Zat(-4)>0} 



ds 



dx > . 



(35) 



Now, since 1 — e ^ < x, for all x > 0, we have that 



l-e-IoMz.t(A)>o}ds ^ / l{z.,(A)>o}ds. 

Jo 
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Therefore, 



1 _ lo ^{Zst(.A)>0} 



ds 



{Zst{A) > 0) ds. 



(36) 



Due to ([55]) and ([55]) , inequality can be written as 



P t 



.-1 



'^{Xs{A)>o}ds > e 1 < (1 - e ^ 



exp I - / / P^ {Z,t {A) > 0) dx ds 

{ Jo JR'i 



(37) 



where 



/ / P^ {Zst{A) >0)dxds= [ [ P^ {Zst{A) > 0) dx ds+ [ [ P^ {Zst{A) > 0) dx ds 

Jo JR'^ Jo J D{st,5) Jo J D{st,SY 

(38) 



with 



Using the inequahty 



D{t, 5) := {x G M'^ : \x\ < 5 > 0. 



[ P^{Zst{A)>0)dx<Kisty 

JD{st,S) 



id/a+y)/2+{l+S)d/a 



which holds for some positive constant K (see Lemma 5 in ^16j), we deduce that 



P^{ZstiA)>0)dxds < K (st)-W"+7)/2+{i+5)rf/"^s 

JD{st,S) Jo 

= ^^-{d/a+7)/2+(l+<5)d/a 



where we used that l — {d/a + j)/2 + {l + 6)d/a > 0. On the other hand, following closely the 
proof of Lemma 5 in [16] one can see that, for sufficiently large t, 

II Px {Zst{A) >0)dxds<ci I P (\\B^\\ > -{st)^/A ds < cit-^. 

Jo JR'i\D{st,s) Jo \ 2 y 

In this way, ([38]) yields the inequality 

"1 







Pa: (ZsM) > 0) dxds < crW»+l)/^+a+S)d/a ^ ^^^-1^ 



(39) 



which is valid for all t large enough, and renders ([55]) . Notice that —{d/a + 7)/2 + (1 + 6)d/a < 
for sufficiently small 6. □ 

Acknowledgement The authors express their gratitude to an anonymous referee for her/his meti- 
culous revision of the paper, and for pointing out a mistake in the proof of an earlier version of 
Theorem [ 
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